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ABSTRACT 

We present the analysis of the bright X-ray binary 4U 1 820-30, based mainly on XMM-Newton-RGS data, but using complementary 
data from XMM-A'ewton-EPIC-pn, INTEGRAL, and Chandra-HETG, to investigate different aspects of the source. The broad band 
continuum is well fitted by a classical combination of black body and Comptonized emission. The continuum shape and the high flux 
of the source {LI Leu ~ 0.16) are consistent with a "high state" of the source. We do not find significant evidence of iron emission at 
energies > 6.4 keV. The soft X-ray spectrum contain a number of absorption features. Here we focus on the cold-mildly ionized gas. 
The neutral gas column density is Nh ~ 1.63 x 10 21 cirT 2 . The detailed study of the oxygen and iron edge reveals that those elements 
are depleted, defined here as the ratio between dust and the total ISM cold phase, by a factor 0.20 ± 0.02 and 0.87 ± 0.14, respectively. 
Using the available dust models, the best fit points to a major contribution of Mg-rich silicates, with metallic iron inclusion. Although 
we find that a large fraction of Fe is in dust form, the fit shows that Fe-rich silicates are disfavored. The measured Mg:Fe ratio 
is 2.0 ± 0.3. Interestingly, this modeling may point to a well studied dust constituent (GEMS), sometimes proposed as a silicate 
constituent in our Galaxy. Oxygen and iron are found to be slightly over- and under-abundant, respectively (1.23 and 0.85 times the 
solar value) along this line of sight. We also report the detection of two absorption lines, tentatively identified as part of an outflow of 
mildly ionized gas (f 0.5) at a velocity of ~ 1200 km s~'. 
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' 1 . Introduction the 10;um feature in ter ms of Mg-rich si licates, if non-spherical 

■ , grain shapes are used dMin et al.ll2007l) . Glassy material, con- 

, The interstellar medium (ISM) in the plane of our Galaxy is a sisting of Mg . rich si licates with metallic iron and sulfide inclu . 

(N| . dynamic and complex environment, composed of mainly neutral sion (called GEMS> | B] . ad i ev | have a i so bee n commonly 

; matter in both gas and dust form and by a warmer gas phase in found during the Stardust mission . Their origin is mostly from 

. the form of diffuse emission in, and above, the Galactic plane. the i nter pianetary environment (e.g. Keller & Messenger 2004), 

The properties of the cold phase in the diffuse ISM have been but a fraction have a composition compatible with an ISM origin 

extensively studied atlong wavelengths, from the far-infrared (Kdler & Messe nger 2008). Iron is however a highly depleted 

to the far-UV (e.g. |DraineJ |2003l for a review). A sizeable element (70 _ 99% of Fe is in dust> | W ilms etaLllS [Whittet 

fraction of the cold phase is locked up in dust grains (e.g. 2 ^ who se inc i us ion int0 sond grains is not completely under- 

b . I Savage & Sembach|1996t | Jenkins | 2009l and references therein). stood ( whittet 2003) Thifj ifj mainly due to the difficulty of 

Amorphous silicate materials together with graphite and poll- modeling iron emiss i n, which does not display any sharp fea- 

cyclic aromatic carbon should account for the majority of the ture m i on g. wave length spectra 

depleted elements measured in the ISM: C, O, Fe, Mg, Si (e.g. ~ ~ 

IWeingartner & Draindl200Tt IWoodenll2008b . One of the major 

spectroscopic signatures of the presence of Fe- and Mg- rich The abundances of the most important metals in the Galactic 
amorphous silicates is the 10/im emission feature. Extensive disk smoothly decrease with the galactocentric distanc e. The av- 
studies of this feature lead to the con clusion that the Fe: Mg ra- erage slope of the distribution is ~ 0.06dexkpc~' (IChen et alj 
tio should be approximatively 1 (e.g. iLi & Draindl200lh . Main [2003, and references therein). However, a large scatter in the 
sources of both Mg and Fe silicates are O-rich asymptotic gi- abundance measurements as a function of the Galactic radius is 
ant branch stars (AGB) and supernovae. However, the process reported. This is attributed to different factors which contami- 
of amorphization of the dust agglomerates, for i nstance by rapid nate the smooth mixing due to the pure stellar evolution pro- 
cooling of the gas phase (iWooden et al.ll2005l) or cosmic rays cess. Indeed the medium can be locally influenced by e.g. su- 
bombardment (ICarrez et al.l 120021) . strongly favors the survival pernovae ejecta, and in-falling m etal-poor material into the disk 
of Mg silicates. In fact, a recent analysis successfully models (Lug aro et al.l lT999: Nittler 2005|). 
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In recent years it has become clear that the X-ray band 
could provide an excellent laboratory to study the silicate 
content of the diffuse ISM, as the absorption K edges of O 
(E = 0.538 keV), Mg (E = 1.30keV) and Si (E = 1.84keV) 
together with the Fe LII and LIII edges (0.71 and 0.72 keV, re- 
spectively) fall in the low-energy X-ray band. The method used 
is to study the absorbed spectrum of bright X-ray binaries, lo- 
cated in different regions in the disk, observed with high-energy 
resolution instruments. This allows to probe the interstellar 
dust (ID) content in a variety of environments, with different 
extinction and with possibly different dust formation history. 
Previous studies of X-ray spectra taken along different lines of 
sight led first to the recognition that not only gas bu t also dust 
plays a role in shaping the iron and oxygen edges (Take i et alj 
2002; iKaastra et all 120091) and later led to the quantitative 
modeling of those edge s dLee et alJ[2009t Ide Vries & Costantinil 
20091: iPinto eTaiM oiO). There is not yet a clear picture of the 
chemical composition of the ID as seen in X-rays. Silicates 
containing andratite (iron-rich silicates ) were reported s tudying 
the oxygen edge of GS 1826-238 dPinto et all 1201 Ol) . while 
iron oxides, rather than i ron silicates we re reported along the 
line of sight of Cyg X-l dLee et aD|2009). This may point to a 
chemically inhomogeneous distribution of ID. 

4U 1820-30 is an extensively studied source, by virtue of its 
extraordinary intrins ic properties. It is an ultracompact (orbital 
period 1 1 .4 minutes, IStella et ai]ll987l) X-ray binary consistin g 
of a neutron star and a He white-dwarf (Rappap ort et al.lfl987l) . 
The presence of X-ray bursts associated with the ne utron star has 
been early recognized (first by Grindla v et al.ll976f). 4U 1820-3 
is classified as an atoll source dHasinger & van der Klisl U989). 
showing kilohertz quasi-periodic o scillation at different frequen- 
cies i n its power spectrum (e.g. Smal e et aDl 1997t IZhang et"aT] 
fl998h . 

The broad band spectrum has been studied with several instru- 
ments. The source displays a Comptonized c ontinuum and a soft 
black body component (e.g. Sidoli et al. 200U). 

With the advent of high resolution spectroscopy, the low- 
energy X-ray spectrum also revealed interesting features. 
Abso rption by ionized ions has been reported in several stud- 



ies (Futamoto et al. 



2004; Yao & Wangl l2005t Uuett et alJ l2006: 



ICackett et al.l 2008a ). The absence of blueshift in the lines of 



this gas and lack of variability genera lly points to an interstel- 
lar origin (but see lCackett et al. 2008a). The hypothesis of a gas 
intrinsic to the source is intriguing as other X-ray binaries of- 
ten display absorption by ionize d gas either at the r est frame of 
the source dvan Peet et alJl2009l) or outflowing (e.g. Mille r et all 
2006] iNeilsen & Ledl2009h . The source has also been used as 
a backlight to illuminate the interstellar dust along the line of 
sight. This results in a dust sc attering halo, which in th is source is 
moderate (~ 3.2% at 1 ke V. iPredehl & Schmitlll 19951) . given the 
relatively low Galactic column density (Nh ~ 1.63 x 10 21 cm" 2 , 
this study). 

Absorption by cold interstellar dust has been never studied 
in detail in this source. In this paper we aim for a comprehensive 
view of the cold absorbing medium along the line of sight to 
this source. This study benefits from the combined information 
provided by the Chandra and XMM-Newton high resolution 
instruments, which allows us to meaningfully study both the 
FeL and OK edges. In addition, we make use of du st models 
updated w ith the latest laboratory measurements (e.g. lLee et alJ 
2008. I2009I) . We also present the broad band continuum behav- 
ior underlying the absorption features. 



The adopted protosolar abund a nces follow the prescrip- 
tion given _by lLodders & Palmel (2009) and discussed in 
Lodders (2010). The broad band energy spectrum (EPIC-pn and 
INTEGRAL) is fitted using the x 1 minimization method, taking 
care that at least 20 counts per bin are present in each data set. 
For the high resolution spectra (RGS and Chandra-HETG) the 
Cash statistic has been usecQ. The errors quoted are for 68% con- 
fidence level, corresponding to A^ 2 = 1 or AC 2 = 1. The spec- 
tral fi tting package used in this paper is SPEX0 ( Kaastra et al. 
1996). The adopted distance is 7.6kpc dKuulkers et alJ l2003). 
The nominal hydr ogen column density toward the source is 
1.29 x 10 21 cirT 2 (Kal berla et alj 120051). to be compared with 
1.52 x 10 21 cirr 2 dDickev & Lockmanlll990l) . 

This paper is organized as follows: in Sect. [2] we illustrate 
the data handling for the different instruments used. Sect. [3] 
is devoted to the continuum determination, using EPIC-pn 
and INTEGRAL data. In Sect. [4] we describe in detail the 
modeling of the different absorption components using RGS 
and Chandra-MEG. The discussion can be found in Sect. [5] and 
the conclusions in Sect. [6] 



2. The data handling 

2.1. Epic-pn 

4U 1820-30 was observed by XMM-Newton for 41 ks on April 
2, 2009 (revolution 1706). The dat a reduction was per formed us- 
ing SAS (ver. 9.0). The EPIC-pn dStriider et alj fcOOl) was oper- 
ated in full-frame masked mode. In this mode, the central 13x13 
(in RAW coordinates) pixels are flagged as bad on board, mask- 
ing the central part of the PSF The full-frame mode was chosen 
to optimally study the halo of diffuse emission surrounding the 
source, caused by the scattering of the intervening dust in the 
line of sight. However, the spectrum of the central source is re- 
covered from the so called Out of Time (OoT) events, which are 
the photons trailed through the detector during the read out of 
the frame. The exposure time of the OoT events is 6.3% of the 
total exposure. We extracted the source spectrum from the RAW 
coordinate event file. We selected 1 1 columns from RAWX 32 to 
RAWX 42, avoiding the central 2 columns, which were affected 
by pile up. We selected RAWY<89 to avoid the innermost region 
of the PSF, which is still affected by pile up or X-ray loadingj. 
The background was selected in a neighboring region. The expo- 
sure time of the background was scaled to match the effective ex- 
posure of the OoT events. Given the brightness of the source, the 
nominal background contribution is modest. However the whole 
area is contaminated by the photons of the scattering halo, which 
have a different distribution depending on the position in the de- 
tector. The OoT events were extracted from a modified event 
file processed in such a way that the charge transfer inefficiency 
(CTI) corrections are not applied and all the photons are consid- 
ered as coming from the central source coordinates. 

2.2. RGS 



Due t o the brightness of the source, some RGS dden Herder et al.l 
1200 lb portions of the grating data were affected by pileup. In 
addition, in RGS2 the high number of events per CCD frame 
(twice the exposure time compared to RGS 1 , due to the single 



1 http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec/manual/XSappendixCash.ht 

2 www.sron.nl/spex 

3 http://xmm2.esac.esa.int/docs/documents/CAL-TN-0050- 1 -0.ps.gz 
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Table 1. The multi-instrument observation log for 4U 1820-30. 



12 14 16 

Wavelength (A) 



Fig. 1. Upper panel: comparison between the 2nd and 1st order 
for RGS1. Lower panel: the ratio between the 2nd and the 1st 
order provides an indication of pile up in the spectrum. 



node readout) tended to fill up the buffers too fast. Limits on 
the data handling in relation to telemetry prevented the buffers 
to be cleared in time, leading to data losses for some CCD's. 
A method to recognize the spectral regions not affected by pile 
up is to compare the first and the second order for each grat- 
ing. Calibration shows that in absence of pile up, the ratio of 
the spectra of the two orders, in the region where they overlap, 
should be unitjQ. In Fig. [TJ we show the ratio for RGS1. We 
see that the maximum pile up is recorded between 14 and 19 
A, where the effective area of RGS is the largest. In this region 
the ratio is about 1.1, which translate in only a 6-7% pile up. In 
RGS2 the ratio between 2nd and 1st order reaches even a value 
of five in selected regions. At longer wavelengths (A > 20 A) the 
combined effect of absorption and decrease of the effective area 
lowers significantly the count rate. Therefore, this region is not 
affected by pile up and can be safely used. Absorption features 
are however in general minimally influenced by this effect. We 
have chosen to keep the useful data given by RGS2 around the 
iron L edge, with care of locally fitting a different continuum 
for RGS 1 and RGS2, using "sectors" in SPEX, as defined below 
(Sect. EJl. 

We also used archival RGS data (Table[TJi in order to improve the 
signal-to-noise (S/N) ratio. This observation displayed a short 
period of flaring background, which was filtered out. This re- 
sulted in a cut of 3 ks on the total exposure time. 

2.3. Chandra-HETG 

In order to gain more signal to noise and energy resolution at 
wavelengths < 18 A, where the Fe, Mg and Si edges stand, we 
considered archival Chandra-HETG data, selecting the higher- 
flux data sets available (Table[TJ. We obtained the final products, 
processed in April 2010, from the TGCat archiv^. The three 
observations displayed similar fluxes and continuum parameters, 
therefore we combined the data, separately for the HEG and the 
MEG arms. HEG data are of low signal-to-noise at wavelengths 
larger than ~17A. In the following we refer to the MEG data 
only, unless otherwise specified. 



Inst. 


orbit/obsid 


date 


net exp. 


rate 1 






(dd/mm/yy) 


(KS) 


(c/s; 


XMM-pn 


1706 


02/04/09 


39.8 


20.9- 


INTEGRAL-JEMX 


790 


03/04/09 


0.61 


28.8 


INTEGRAL-1SGRI 


789 


31/03/09 


42.4 


15.5 


XMM-RGS 


1706 


02/04/09 


41.5 


27.6 


XMM-RGS 


0336 


09/10/01 


34.7 


33.1 


Chandra-MEG 


6633 


13/08/06 


25 


125.1 


Chandra-MEG 


6634 


20/10/06 


25 


172.8 


Chandra-MEG 


7032 


05/11/06 


46 


146.2 



4 http://xmm.vilspa.esa.es/docs/documents/CAL-TN-0075- 1 -0.pdf 

5 http://tgcat.mit.edu/ 



Rates refers to the full band for each instrument. 
2 Rate from selected pixels of the OoT events. 



2.4. INTEGRAL 

In order to constrain the continuum spectrum of 4U 1820-30 
across an as-large-as-possible energy band, INTEGRAL data 
can be very useful. Fortunately, the Galactic center region was 
observed by INTEGRAL for ~ 73 ks from March 30, 2009 to 
March 31, 2009 during satellite revolution 789, and for ~ 38 
ks on April 3, 2009 during revolution 790. The XMM observa- 
tion ended during the latter observation. Therefore we can con- 
sider these INTEGRAL observations as quasi-simultaneous with 
the XMM -observation. 4U 1820-30 was in the field of view of 
the INTEGRAL So ft Gamma-Ray Imager ISGRI (15-300 keV; 
Lebru n et al . 2003) at source angles smaller than 14°5 for all 21 
individual pointings (so called science windows) of the rev. 789 
observation and all 10 pointings of rev-790. The source was out- 
side the field of view (13°2 diameter zero response) o f the Joint 
Europ ean Monitor for X-rays JEM-X (3-35 keV; iLund et al.1 
120031) during the rev. 789 observation, but within the field of 
view during 4 individual pointings of rev. 790. 

We generated mosaic maps for 10 logarithmically spaced en- 
ergy bands across the 20-300 keV range for the ISGRI rev. 789 
observation using OSA version 9.0 (distributed by the ISDC; 
ICourvoisier et al.ll2003l) imaging tools. We checked these maps 
for significantly detected sources and subsequently derived spec- 
tral information for all these sources in 13 pre-defined energy 
windows across 13-520.9 keV. 

For JEM-X (telescope 1) we followed similar imaging and 
spectral extraction procedures yielding spectra of 4U 1820-30 in 
16 pre-defined energy bands across the 3.04-34.88 keVThe ob- 
servation log is displayed in Table [TJ 



3. The broad band spectrum 

We first extracted the light curve from the EPIC-pn OoT events 
as described in Sect. 12. II therefore the resulting count rate cannot 
be used to recover the absolute flux. The 0.5-10keV light curve 
is shown in Fig. [2] using a bin size of 500 s. It shows a modest 
variation of at most 13% during our observation. The rates in 
the soft (0.5-2keV) and the hard (2-10keV) band followed the 
same pattern. 

We fitted the EPIC-pn and INTEGRAL spectra simultaneously. 
The source was caught at a relatively high flux, with a 2-10keV 
flux of ~ 5.0x 10~ 9 erg cirT 2 s _I (Tabled. In Fig. [3](upper panel) 
we show the residuals with respect to a simple powerlaw with 
photoelectric absorption. From the residuals we see a complex 
spectrum at energies below 1 keV, probably due to ionized ab- 
sorbing gas. The spectrum is well detected up to 40keV. 
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TIME (s) 
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Fig. 2. Light curve of EPIC-pn data as selected in Sect. I2.ll Note 
that the displayed count rate cannot be used to recover the flux 
of the source. The bin time is 500 s. 

Table 2. Broad band modeling of the source using EPIC-pn and 
INTEGRAL data. 



n: 



tst - 

H 
Tbb 

kT 

x 2 /y 

F 0.5-2 
-F2-10 



7.8 ±0.3 
0.26 ± 0.01 
0.48 ± 0.02 

3.9 ±0.1 
5.4 ±0.2 
1756/1720 

1.32 ± 0.02 x 10- 9 
5.02 ± 0.07 x 10- 9 
2.70 ± 0.02 x 10" 9 



Notes: 
Units are: 



10 cm - for column densities, keV for temperatures, Fluxes 



are measured in erg cm s . 

1 from RGS data 

2 from JEMX data 

3 from JEMX and ISGRI data. 



As a simple powerlaw is not an acceptable model (^ 2 /v = 
2412/1724 = 1.39, where v is the number of degrees of free- 
dom), we added first a black body component, to mimic the often 
observed soft energy emission (x 2 /v = 2273/1722 = 1.32). This 
is parameterized by the black-body temperature (I'm) and its 
normalization. However a black body plus power law model fails 
to explain the spectral curvature seen at high energies. We then 
substituted the powerla w model with a C omptonization model 
(COMT model in SPEX: lTitarchukl 19941) . where the parameters 
are: the temperature of the seed photons (kTu), the optical depth 
of the electron cloud (r) where the photons are Compton scat- 
tered to the final temperature (kT\). This model provides a sat- 
isfactory fit to the broad band continuum Of 2 /v = 1756/1720 = 
1.02; see Fig. [3] lower panel, Fig. [4] and Table|2|. 

We repeated the fits using a disk black body model (i.e., 
emission from a standard Shakura-Sunyaev disk, model DBB 
in SPEX) instead of a simple black body. We do not obtain 
an equally good fit using DBB either in addition to a simple 
powerlaw model (x 2 /v = 2302/1722 = 1.33) or DBB plus 
Comptonization ( X 2 lv = 1998/1720 = 1.16). 

In the 6.4 keV region, there is no clear evidence of iron 
emission lines. A 3 <x upper limit of 25 eV on the equivalent 




+ 



i^mmm^ ^h- 



Observed Energy (keV) 

Fig. 3. Upper panel: residuals to a simple absorbed power law 
for EPIC-pn and INTEGRAL data. JEMX data have been high- 
lighted for clarity. Lower panel: residuals to the best fit con- 
sisting of a black body plus Comptonized emission, modified 
by both a neutral and ionized absorbing material. EPIC-pn data 
have been rebinned for displaying purpose. 




Observed energy (keV) 



Fig. 4. Upper panel: The best fit model consists of a black 
body (dashed line) plus Comptonized emission (dot-dashed 
line), modified by both a neutral and ionized absorbing mate- 
rial (Table 0. Lower panel: Best fit model superimposed to the 
EPIC-pn, JEMX and ISGRI spectra. 



width is obtained if a delta line is put at a fixed energy of 6.4, 
6.70 and 6.97 keV, respectively. We also searched for a rela- 
tivisti cally modified line profile (LAOR model in SPEX; lLaorl 
11991b . In this model, the line arises from the accretion disk 
where the intensity follows a R q profile, where R is the dis- 
tance of the disk gas from the source. Leaving the energy as 
a free parameter did not lead to a meaningful fit. We then fix 
the energy of the line at 6.97 keV, fol lowing previous studies of 
4U 1820-30 (e.g. lCackettetalJl2010h . The best fit yelds a very 
broad, but weak line profile, with parameters q = 2.8 + 0.2 and 
i < 14°, where i is the disk inclination. The flux of the line is 
(1.6±0.7)x 10~ 4 photons cm^s" 1 . 
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The complex absorption at soft energies has been mod- 
eled following the prescription given by the high-resolution data 
(Sect 12.21 ). Therefore we added a minor contribution from ion- 
ized absorbers (A^ n ~ few x 10 20 cm" 2 ) and the Galactic neu- 
tral absorber, which contributes to most of the low-energy cur- 
vature. We note that the best fit neutral column density is about 
40% lower than the best fit obtained in the RGS analysis. This 
may be caused by foreground scattering halo soft-emission (ob- 
served for the full exposure) located just in front of the OoT 
event (observed for 6.3% of the exposure time) which still re- 
mains after the subtraction of the local background. The scat- 
tering halo appears indeed as diffuse emission which extends 
on top of both the wings of the point spre ad function of the 
source and the OoT events themselves (e.g. iPredehl & Sc hmitt 
119951: ICostantini et al.ll2005[). The scat tering process is strongly 
energy dependent dMathis & Ledll99ll) . In particular, the halo is 
brighter at softer energies. The net effect is therefore to add more 
photons to the soft energy spectrum of the OoT events, reducing 
the measured absorption toward the source. This effect appears 
even more enhanced by the reduction of the OoT absolute flux 
due to the cut of the central two columns, where most of the 
source photons are. 

4. The high resolution spectrum 

Since the energy band of the RGS is insufficient to discriminate 
among different broad-band models, we adopted the model de- 
fined in the previous section. The spectrum shows evident O i K 
and Fe i L edges, due to absorption by neutral material. In addi- 
tion, absorption by ionized gas is highlighted by the O vm, O vn 
and O vi absorption lines. 

We used a collisionally ionized plasma model (model HOT in 
SPEX), with a temperature frozen to kT = 5 x 10~ 4 keV, in order 
to mimic a neutral gas. This model fits to first order the edges 
and resonant lines from the neutral species, leaving however no- 
ticeable residuals in the fit. In order to extract more robust results 
from the absorbtion features, we added archival RGS (Sect. \2.2\ 
and Chandra-MEG data (Sect. [273b . As these data sets were ob- 
tained at different epochs, the continuum shape may differ sig- 
nificantly. To bypass this complication, we used the "sectors" 
option is SPEXO and we left the continuum parameters free to 
vary for each sector. Note that the RGS1 and RGS2 of the re- 
cent observations were treated as different data sets (thus with 
different continua), as the RGS2 broad band shape was affected 
by pileup. The parameters of absorber components were coupled 
together for all the data sets. 

4.1. The ionized gas 

The RGS spectrum of 4U 1820-30 displays narrow absorption 
features from ionized gas, mainly from oxygen, iron and neon. 
In particular Ne ix, O vi-O vm are prominent features of this gas. 
For this fit we ignored the data below 13 A for RGS. In this re- 
gion some wiggles in the local continuum, mainly caused by pile 
up, could add uncertainties in the determination of the absorp- 
tion parameters. Here we focus on the neutral or mildly ionized 
phase of the ISM. Therefore we simply model the higher ion- 
ization lines with a phenomenological model (SLAB in SPEX) 
which calculates the transmission from a thin layer of gas. We 
defer a more detailed analysis of this gas component to a future 
publication. The resulting ionic column densities of the main 

6 See Ch. 5 of the SPEX cook book: 
http://www.sron.nl/files/HEA/SPEX/manuals/spex-cookbook.pdf 



Table 3. Parameters of the main lines of the more ionized ions 
as measured by RGS. 



ion 






cm 2 




Cvi 


16.7 ± 


0.1 


Nvi 


< 15.1 




N vn 


15.5 ± 


0.4 


Ovi 


15.4 ± 


0.3 


Ovn 


16.2 ± 


0.1 


vm 


17.6 ± 


0.1 



Note: 

The assumed velocity broadening is <r r = 62 km s 



ions, as measured b y RGS, are reported in Table [3] Following 
lYao & Wa ng (2005), we assume a velocity broadening of cr v = 
62 km s _1 . All subsequent fits in the present analysis already take 
into account this highly ionized component. 

4.2. The oxygen edge 

The oxygen edge at 0.538 keV has been fitted using the high 
quality data of the two RGS epochs between 19-36 A. The neu- 
tral material, modeled by a low temperature (kT — 5 x 10 -4 keV) 
in a collisional ionization equilibrium, fits well the long wave- 
length curvature due to the ISM absorption as well as the Oi 
ls-2p absorption line at 23.5 A (Fig.|7]i. The edge shape appears 
modified by several sharp absorption features. Some of them 
can be easily identified with the O vn and O vi lines, belonging 
to a highly ionized phase (Sect. 14. 11 1. At the shorter wavelength 
side of the Oi line, we see On at A — 23.35 A and also weaker 
Om and Oiv lines. Two absorbers with different temperatures 
are required to fit these lines (Table|4] column (1 )). The reported 
line widths are automatically evaluated in SPEX using a curve of 
growth analysis applied to multiple l ines belonging to the same 
absorption system (e.g. Spitzerl [l978l) . This allows to accurately 
evaluate the line width using basic parameters such as the ionic 
column density and equivalent width. In this case, Oi, Fei, 
Ni transitions are used for the cold gas (comp 1 in Table |4j. 
The strongest transitions are saturated, which in principle may 
introduce additional uncertainty in determining the line width. 
However, as shown in Fig. [5] for the lines we study here, the 
ionic column density as a function of the line equivalent width 
depends marginally on the velocity width. This behaviour of 
the curve of growth is due t o the large a - Voigt parameter for 
the inner-shell transitions (Mihalas 1978: iKaastraet al.ll2008l 
for a full discussion). We note that the nitrogen region is not 
modified by dust and the line widths are then easier to evaluate. 
In addition, many of the the weaker transitions of the same 
neutral ions are not saturated. For the kT ~ 3.2 eV gas (comp 2), 
O n, Fe ii, N ii are among the strongest ions. These are relatively 
weaker lines where saturation does not play a major role. The 
relatively limited resolution of the instrument does not allow 
to distinguish multiple velocity-width components within a 
same absorption line. Therefore we consider here the total ionic 
column density along this line of sight. 

Dust should also be present along this line of sight, 
as 4U 1820-30 is know n to display a dust scattering halo 
(IPredehl & Sch mitt 1995), which is also evident in the present 
EPIC-pn data. We therefore fitted the oxygen edge adding the 
AMOL model in SPEX. We considered 16 oxygen compounds 



5 



E. Costantini et al.: XMM-Newton observation of 4U 1820-30 




Fig. 5. Line equivalent width as a function of the ionic column 
density for O i, N i and O n. The velocity dispersions are <x,, =10, 
100, 200km s . Curves of different ions have been shifted for 
clarity. 



for this analysis (see Appendix A). This list includes silicates 
and oxides which are common in the ISM plus some lighter 
(also icy) materials. Absorption features by dust, contrary to 
gas, are smooth and broadened. All the oxygen compounds con- 
sidered here have no appreciable effect at A £ 23.7 A (Fig. IA.ll ). 
Below we describe a simple fit using a single dust component 
at a time, adding a second component if the fit requires it. This 
makes the description of the model and of the steps to the final 
best fit simpler to follow. This approach is however the result of 
a rigorous fitting procedure for both the oxygen and iron edges, 
which is fully described below (Sect. l4~4l >. 

We first attempted to fit the whole spectral region only in 
terms of dust absorption. In the fit we let the hydrogen column 
densities of the gas producing Oi and On (comp 1 and 2 in 
Table |4]i as free parameters. We also let the abundance of oxy- 
gen in the cold phase (comp 1) as a free parameter, in order to 
balance the amount of oxygen locked up in dust. In this model, 
we suppressed the gas component (comp 3) mainly producing 
Om (A = 23.07 A), as this line is the main absorption structure 
in the region where also dust plays a major role. It is clear that 
fitting the region around 23 A only in terms of dust remarkably 
worsens the fit. The worsening range is indeed AC 2 = 35 - 180, 
depending on the compound, for Av = 2, with respect to the 
pure gas modeling. This indicates that not only dust is absorb- 
ing the spectrum in the ~ 22.7 - 23.1 A region. In Fig.|6l upper 
panel, we show the oxygen region where the pure gas fit is com- 
pared with the fit with MgSiO^ (dashed line) or Mgi eFecuSiCU 
(olivine, dashed-dotted line), suppressing comp. 3. 

Finally, we included both the kT ~ 13eV gas absorption 
(comp 3) and dust absorption (Table. [4] column (2)). The 
improvement of the fit is then AC 2 = 10 - 108 for Av = 1. Most 
of the individual dust components improve the fit. However, 
we can identify two compounds for which the improvement of 
the fit is the largest (AC 2 = 105 - 108). These are H 2 0, in the 
form of crystal or amorphous ice, and MgSiC>3. In Fig.[6](lower 
panel) we show how a fit with pure gas and a gas+MgSiC>3 
mixture (comp 5 in Table |4]i compare with each other. Note that 
a fit using H2O is non distinguishable in practice from MgSiOi 



Fig. 6. Fits in the oxygen region. Here only RGS 1 of the most 
recent observation is displayed for clarity. Upper panel: com- 
parison between a fit only in terms of gas (solid line) and the 
best fit obtained using dust (dashed line). For comparison also an 
unacceptable fit using olivine (dashed-dotted line) is displayed. 
Lower panel: comparison between a fit only in terms of gas 
(solid line) and the best fit obtained with a mixture of gas and 
dust (dashed line). See more details in the text. 

as the spectral feature is very similar (Fig. IA.1I Sect. 14.41 for 
discussion). Interestingly we note that the compounds with 
both oxygen and iron provided the least improvement (if 
not a worsening) of the fit. In particular, the fit rules out the 
most complex aggregates of oxygen and iron (e.g. magnetite, 
franklinite, olivine and almandine, see Fig. lA.ll for the chemical 
composition). 

As noted above, none of the dust compounds considered 
here have features below 23.7 A. However, from Fig. |6]we note 
two clear absorption features (present in both RGS data sets) at 
~ 22.3 and 22.65A (marked with arrows) detected with 6.8 and 
3.8<t significance, respectively. These features do not belong to 
the ionized absorber producing the O vn line (Sect. l4~TI ), but to a 
lower ionization gas. We tentatively fitted these features with a 
photoionization model (XABS in XSPEC), obtaining a low col- 
umn density (A^h ~ 4.3 x 10 19 cirT 2 ) and a ionization parametefl 
log£ ~ -0.5 (comp 6 in Table @J. Interestingly, the two major 
lines, identified as O iv and O v present a systematic blue-shift 
of about -1200km s _1 (Table|4] column (3)). From the absorber 
model we notice that the other lines predicted using these pa- 
rameters fall either in a lower-resolution and noisier part of the 
spectrum (e.g. lines of iron) or where the spectrum is heavily ab- 
sorbed by the neutral gas (e.g. nitrogen lines). 
The best fit in the oxygen edge region, including RGS data sets 
taken in two epochs is displayed in Fig. [7] 

4.3. The iron L edge 

In Fig. [8] a fit of the Fe LII and LIII edges in terms of absorption 
by pure gas with solar abundances is shown (dotted line). In this 
fit two RGS epochs and the MEG data sets were used. There is 

7 The ionization parameter £ is defined as £ = L/nr~, where L is the 
ionizing luminosity, n the gas density and r the distance of the gas from 
the source. 
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Table 5. Goodness of fit for the iron edge models. 



Fig. 7. Best fit transmission spectrum (Table [4] column (3)) in 
the oxygen edge region. This allows to compare absorption at 
different epochs, removing the contribution of the continuum. 
The upper curves display the transmission of the various absorb- 
ing components. Note that for clarity only absorption by oxygen 
is displayed (i.e. no hydrogen absorption). Data have been re- 
binned for clarity. 

Table 4. Parameters for the oxygen edge modeling. 



model 


C> 


Gas (A Fc fix. ) 


2377/981 


Gas (Ap c free) 


1544/980 


Fe 2 Si0 4 


1505/977 


FeS0 4 


1486/977 


FeO(OH) 


1464/977 


Fe 2 3 


1442/977 


Fe met. 


1424/977 


Fe met.+Fe203 


1421/976 







(1) 


(2) 


(3) 


comp 1 




19.6 ±0.2 


17.3 ±0.2 


16.3 ±0.2 




jcold 


0.5 fix. 


0.5 fix. 


0.5 fix. 






20 ± 10 


20 fix. 


20 fix. 


comp 2 


N H 


0.50 ± 0.09 


0.56 ± 0.09 


0.61 ±0.09 




O-y 


< 10 


10 fix. 


10 fix. 




T 


3.2 ± 0.4 


3.2 ± 0.4 


3.3 ± 0.4 


comp 3 


N H 


0.90 ± 0.08 


0.36 ± 0.06 


0.26 ± 0.06 




T 


9±2 


13 ±3 


13 ±3 






100 fix. 


100 fix. 


100 fix. 


comp 5 


N 




2.2 ± 0.2 x 10- 3 


2.3 ± 0.2 x 10~ 3 




N si 




7.4 ± 0.8 x 10- 4 


7.8 ± 0.8 x 10- 4 




N Ms , 




7.4 ± 0.8 x 10~ 4 


7.8 ± 0.8 x 10~ 4 


comp 6 


N s 






0.43 ± 0.06 




logf 






-0.5 ±0.1 




"out 






-1200 ± 130 




C 2 /v 


5863/4572 


5745/4571 


5684/4568 



Notes: 

Units are 10 20 cirr 2 for column densities (N H ), eV for temperatures (T) 
and kms~' for line broadening (<r v ) and the outflow velocity (v ou t)- 
Column(i): model with neutral and mildly ionized gas. Column (2): 
mildly ionized gas + dust. Column (3): milidly ionized gas + dust + 
mildly ionized outflowing gas. In all models the highly ionized compo- 
nent is already included. 



a clear mismatch not only in the depth, but also in the position 
of the of the LII and LIII edges. We note that in the fit of the 
iron edges, it is essential to include the higher resolution MEG 
data (Fig. [8]i. A straightforward way to improve the fit is to mod- 
ify the iron abundance. Letting this parameter free, we obtain 
a reasonable fit, although the position of the edge in the model 
still does not match the data. In this case the abundance of iron, 
assuming that the absorption is only in the gas phase, is about 
0.37 times the solar one. However, absorption by dust is known 



to alter the LII/LIII ratio and to shift the position of the edge 
with respect to absorption by gas (see e.g. lLeeetal. I l2009h. For 
the fitt ing, we considered the compounds measured by Lee et al.l 
(2009): metallic Fe, hematite (Fe 2 3 ), lepidocrocite (FeO(OH)), 
fayalite (Fe2Si04) and iron sulfate (FeSOzt). These transmission 
models have been recently implemented in the AMOL model in 
SPEX. As for the modeling of the oxygen edge, we first followed 
a rigorous approach (Sect. 14.4b which in turns justifies the sim- 
pler approach described here, i.e. fitting the dust components one 
by one, together with the gas model. In the latter, the iron abun- 
dance is a free parameter. This procedure provides a measure 
of the depletion of the gas phase. On the goodness of fit bases, 
the single dust compound which best models the data is metallic 
iron (Table [5}. Then we tried all possible combinations of gas 
plus two dust components. We only find a marginally significant 
presence of Fe203, in addition to the metallic iron component 
(Table [5] Fig. [8] solid line). Hematite (Fe203) is a compound 
for which both iron and oxygen edges' dust profiles are imple- 
mented in the AMOL model. We therefore tested the presence 
of this compound in the oxygen edge region, obtaining an upper 
limit on the column density, consistent with the value obtained 
in the iron region. In Fig. [8] we also show the comparison be- 
tween the best fit model and a model with an iron rich olivine 
(fayalite, Fe2Si04), which does not provide a good fit (Table[5]l. 
This suggests that fayalite cannot have a major contribution to 
the absorption here. We cannot test the same kind of olivine as 
for the oxygen edge (Mgi 6Feo.4Si04) as at present no laboratory 
measurement of this compound at the iron L edge is available. 
In Fig. [9] the best fit is shown where data from both RGS and 
Chandra-MEG are displayed. The spectra have been normal- 
ized to their continuum shape for displaying purposes. We also 
normalized to the total hydrogen column density so that only 
absorption by iron is visible for the neutral phase. Along with 
the best fit (solid line), also the different contribution are shown: 
atomic gas (dotted line) and metallic iron and Fe203 (light dash- 
dotted lines). Also the contribution of the ionized gas to the spec- 
trum is shown with a dashed line (e.g. O vn line at 17.76 A). 
Despite the significant improvement of the fit with respect to 
a gas-only model, the iron edge is clearly more complex that 
our parameterization. In particular, positive residuals are present 
on the longer wavelength side of the edge. T his effect has 
been found previous ly using different instruments (Kaastra et al.l 
2009l lLeeetaTl200^) . 



4.4. Further insight in the FeL and OK edge fitting 

The modeling of photoelectric edges modified by dust absorp- 
tion requires some caution, as in principle many components can 
contribute to the edge shape, which is in turn smeared, contrary 
to the atomic transitions features, which are very sharp and 
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Fe 2 Si0 4 + GAS 
• Fe + Fe 2 3 + GAS 



17.0 17.5 18.0 18.5 

Wavelength (A) 



Fig. 8. Detail of the iron L edge region. Here for clarity we dis- 
play only the MEG and a displaced RGS2 data set. The best fit 
(solid line, see Table |5]l is compared with pure gas fit with solar 
abundances (dotted line) and with a mixture of gas and iron-rich 
olivine (thick solid light line). 




16.8 17.0 17.2 17.4 17.6 17.8 18.0 
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Fig. 9. Best fit transmission spectrum of the region around the 
iron LII and LIII edges. This allows to compare absorption at dif- 
ferent epochs, removing the contribution of the continuum. The 
upper curves display the transmission of the various absorbing 
components. 

relatively easy to identify. Therefore the accuracy of the result 
depends on the resolution of the instruments, the signal-to-noise 
ratio, the completeness of the dust models and finally on how 
many independent parameters there are in the problem. 
In this study we combined the higher resolution of Chandra- 
MEG with the sensitivity of RGS to study the Fe L edges. For 
the oxygen edge we used the high quality RGS data, as the 
Chandra-HEG effective area dramatically drops in that spectral 
region. 4U 1820-30 is also a high flux source, absorbed by a gas 
and dust column density which is optimal in order to study the 
O and Fe edges. 

The dust edge profiles that we use in our study are a collection 
of what is available in the literature up to now ( e.g. Barrus etaLl 
1 19971 iParent et al. Il2002t IvanAken et aDll998b iLee et al.l i2008. 



120091) . This data-base is not complete, but the main representa- 
tives of oxides and silicates for both oxygen and iron are present 
(see Appendix A, ILee et al.ll2.Q09b . This is sufficient for a first 
order analysis. 

The simultaneous analysis of two edges produced by dust on 
the same line of sight allows to cross-validate our results. In 
addition, common knowledge on abundances and physical 
conditions in the ISM allow us to put further constraints to the 
edges fitting. Here we illustrate this in detail. 

In the ISM a mixture of chemical compounds exists. We 
tried to simulate this by mixing all possible combinations of 
edges profiles available for a given edge. The AMOL model 
allows to fit four different components at a time. We therefore 
had for each edge a number of fits n where n = nd\j\\(nd - 4)!. 
Here nd is the number of available edge's profiles. Examining 
the output we notice that the combinations which improve the 
fits prefer one or two components, putting the remaining column 
densities to zero. If two components are preferred, one strongly 
dominates over the other, whose relative value is never more 
than 20%. This applies to both the oxygen and the iron edge. 
This justifies the simplified way of illustrating the fitting used in 
Sectg3]|43]above. 

We further analyzed all fits that were within 3 sigma the best 
fit values shown in Sects. I4~2l and l4. 31 Interestingly, for both iron 
and oxygen, no models containing iron-rich silicates are among 
these selected models. 

For iron we could rule out models, within the "3 sigma" set, 
containing hematite (Fe203) with a total molecular column den- 
sit^ larger than 4.4 x 10 15 crrT 2 . This resulted from the cross- 
validation we performed with the oxygen edge, as for hematite 
we have both the O and Fe dust profiles. Indeed, none of the 
"3 sigma" models for oxygen contained this compound. We 
therefore added it to the oxygen edge best fit model, obtaining 
the limit above (and compatible with what found for the iron 
edge best fit, see Sect. 14.31 ). 

We could also rule out the "3 sigma" models containing Fe2S2 
as for this case the fit delivered an amount of sulfur ~3.8 times 
solar. The sulfur edge would have been significantly enhanced, 
but this is not evident in our data. 

For the oxygen edge this multicomponent fitting pointed very 
clearly to a predominance of MgSi03 in all the "3 sigma" set. 
As noted in Sect. 14.21 a similarly good fit is obtained if water ice 
is the predominant component. This is because of the very sim- 
ilar edge profiles of the two compounds (Fig. IA.lt . However, 
we could readily rule out this option, as along the line of sight 
of 4U 1820-30 the bulk of cold material is diffuse interstellar 
medium, where water i s virtually abse nt and silicates are on the 
contrary abundant (e.g. Wooden 20081). 

4.5. The magnesium and silicon edges 

The Mg K and Si K edge (at 9.47 and 6.723 A, respectively) 
were analyzed using the combined information of Chandra- 
HEG and MEG, which provide in this band a higher energy 
resolution. In the case of Si, the energy of the edge is also out- 
side the RGS range. The quality of the data and the relatively 
low Galactic column density hampered any detailed characteri- 
zation of the edges. We used the auxiliary information provided 
by the oxygen edge fitting (Sect. 14.2b where a model includ- 

8 Note that this translates in N < 1.3 x 10 I6 cnr 2 and N < 8.8 x 
10 15 cirT 2 for O and Fe respectively. 
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ing MgSiC>3 was preferred. For that compound we could ben- 
efit of the laboratory measurements of both the O and Si edge. 
For Mg, for which we did not have any measurement, the rel- 
ative column density was included in the model in the form of 
an artificial gas-shaped edge. We therefore added to the gas ab- 
sorption model this specific dust component, fixing the amount 
of magnesium and silicon locked in dust and letting the amount 
of Mg and Si in the gas phase as free parameters. This model 
agrees with the data, within the uncertainties. The column den- 
sity of Si in dust is A^ ust ~ 7.8 x 10 16 cirT 2 , while the Si in 
gas is iV|? s < 1.3 x 10 16 cirT 2 . Therefore the inclusion of Si 
in dust is ~85% of the total neutral Si along the line of sight. 
For the Mg edge, we obtain a N^l < 3 x 10 15 cirT 2 , while the 
value imposed by the oxygen fit to the Mg locked in dust is 
= (7-8 + 0.8) x 10 16 cirT 2 . The amount of Mg locked up in 
dust according to this model is about 96%. 

4.6. The neon and nitrogen edges 

In order to study the neon edge (A = 14.24 A), we used the com- 
bined data of RGS and Chandra, while for the nitrogen edge 
we used only the RGS data. These elements ar e only present 
in gas form in the ISM (e.g. IWilms et al.ll2.Q00l) . Fitting neon 
with our Galactic absorption model we find a column density of 
2.07 + 0.02 x 10 17 cm -2 , which converts in a slight formal over- 
abundance with respect to solar 1.2 + 0.1. 
Next to the nitrogen edge at 30.76 A, we can also distinguish 
absorption by the ls-2p line at 31.28 A. The measured column 
density of nitrogen is (1.33 +0.01) x 10 17 crrT 2 , corresponding to 
an abundance of 0.9 + 0.1 times solar, consistent with the value 
expected using solar abundances. 



5. Discussion 

5. 1 . The broad band continuum and the iron emission line 

We have analyzed the simultaneous observation of 4U 1820-30 
using XMM-EPIC-pn and INTEGRAL. The fit of the continuum 
is consistent with black body emission plus a Comptonized com- 
ponent, the latter extending up to 40keV. The source was ob- 
served in a high flux state, accreting at Lboi/^Edd ~ 0.16. We 
have assumed Lboi = 2 x L(2-20 keV> where the unabsorbed lu- 
minosity Lp-20 keV) is ~ 3.1 x 10 37 ergs~'. The Eddington lu- 
minosity f or ultracompact syste ms is LEdd ~ 3.8 x 10 38 ergs _1 , 
following Kuu lkers et al.l d2003l) . The spectral shape is also con- 
sistent with a high state scenario. The hard energy portion of the 
spectrum is well fitted by a Comptonization model with an ef- 
fective cut-off at ~ 40keV, consist ent with previous observations 
of the source in the high-state (IB loser et al J 120001: [Sidoli et al. 
120011) . The black body emission, which fits well the soft energy 
range of the present observation, has a temperature similar to 
that of the comptonization seed photons. This may suggest that 
in ultra-compa ct X-ray binaries s eed photons originate from the 
accretion disk (Sidoli et al.ll2001l) . 

A relatively weak or undetected iron line is a common fea- 
ture o f the confirmed ultra -compact X-ray binaries (as defined 
e.g. in lin't Zand et al.ll2007l) . as 4U 1820-30 is. This maybe sim- 
ply explained by the nature of the donor star in such systems. 
Due to the high surface gravity (log g ~ 8), heavier elements 
in the wh ite-dwarf companion sink in the internal layer of the 
star (e.g. lFontaine & Michaudlll979l) . leaving little high Z met- 
als for the accretion flow. In 4U 1820-30 however, a faint, vari- 
able and relativistically smeared iron line as been reported (e.g. 



ICackett et al] l2008bl 12010). In our data, no narrow emission 
features, either from neutral or ionized iron, have been found 
(Sect. 0. A fit with a relativistically smeared iron line leads to 
a very elusive 2-sigma detection of the feature. The photon flux 
of this line is about 16 times fa inter than that found in recent 
Suzaku data dCackettet al.ll2010l) . 

In exceptional circumstances, like superburst episodes, a large 
quantity of heavy elements ashes may be ejected during the neu- 
tron star photosphere expansion. This could possibly cause de- 
tectable features (e.g. in absorption) by e.g. iron in the X-ray 
spectra din't Zand & Weinberg! 12010) . In 4U 1820-30 the pres- 
ence of the iron line could then be the tracer of burst activity, not 
detectable at the epoch of our observation. 



5.2. The dust components 

The high resolution spectrum of 4U 1820-30 shows evidence of 
absorption by dust, especially around the Fe L and O K edges 
(Figs. [7] |5). In the fits we explored all the possible combina- 
tions of dust mixtures (Sects. 14.31 and l4~2l . The best fit shows 
a preference for absorption by enstatite (MgSi03), which is the 
end series of the pyroxene silicate mineral, metallic iron (Fe) 
and traces of iron oxides in the form of hematite (Fe203). In 
the modeling, 16 and 5 different compounds were considered 
for oxygen and iron, respectively. The set of models does not 
cover the whole range of possible compounds present in the 
ISM. However what are commonly believed to be the main con- 
stituents (e.g. olivines, pyroxene, oxides as well as simpler com- 
pounds) in the diffuse ISM (e.g. Whit tet1l2003t IWoodenll2008l) 



Lee et alJ2009l) . A further ef- 



are represented (see Appendix A, 
fort in obtaining new measurements has recently been carried out 
(iLee & Rave 112005c iLee et alJ2009L and de Vries et al., in prep.). 
The available X-ray measurements were mostly performed on 
crystalline materials. The presence of dust in the form of crys- 
tals is significantly present only in specific astronomical environ- 
ments, like com ets (|Woodenl2008]) or th e inner regions of proto- 
planetary disks dvan Boekel et alj|2004l) . In the ISM the amount 
of crystalline silica tes, compared to amorphous grains, is < 5% 
dLi & Draindl200ll) . therefore we do not expect to find a size- 
able amount of it in our data. However, for the handful of com- 
pounds for which X-ray laboratory measurements of both glassy 
and crystalline forms are available, we tested that the spectral 
difference is only appreciable if the resolution is of the order 
of 3 - 5 eV. Therefore, with the present data, our analysis using 
crystalline grains is a good first order approximation of model- 
ing the chemical composition of ID. In the glassy form, crys- 
talline material loses the ordinate internal structure proper of a 
crystal. For simplicity, although formally not correct, we may 
call an amorphous silicate with a certain (for instance enstatite) 
stoichiometry, glassy-enstatite or glassy-MgSi03. 



5.3. Abundances and depletion 

As X-ray spectra do not display any sharp H feature, either in 
emission or absorption, the total hydrogen column density has 
been evaluated from the low energy curvature of the spectrum, 
which extends down to 36 A. This is a reliable method, as start- 
ing already at about 25 A, the transmission is largely dominated 
by He and H. The best fit value for the hydrogen column density 
is N H = (1-63 ± 0.02) x 10 21 cirT 2 , slig htly larger than the 
nominal values (1.52 + 0.07) x 10 21 cm" 2 dDickev & Lock man 
1990) and (1.32 + 0.05) x 10 21 cm 2 (Kalbe rla et al.l 120051) . 
which are the average values over a region of 1 radius around 



9 



E. Costantini et al.: XMM-Newton observation of 4U 1820-30 



the source. The reason could be that the spectral curvature 
measures the total hydrogen (i.e. H i, H n and H2), which should 
be larger than H 1 alone . The amount of H 2 in the diffuse ISM 
is relatively low (e.g. iTakei et aD l2002i) . We estimated the 
amount of Hn using the On/Oi ratio a s a proxy (following 
the expression in Field & Steigmanl ll971l) . We derive that ~4% 
of the measured hydrogen column density is in the form of 
Hn. This would reconcile our value with the Hi measurement 
mentioned above. Besides, it could be that a small fraction (in 
this case no more than 1 - 3 x 10 20 cirT 2 ) of the total H resides in 
the imme diate surroundings of the source, like in o ther sources 
(e.g. iPredehl & Schmitd[l99"l Ivan Peet et al.ll2009l) . Such a low 
column density would be insufficient, however, to produce deep 
absorption features. We compute the following abundances 
taking as a reference the value we measure from the X-ray 
spectrum. 

In the iron edge region (Sect. l4~3l >. a fit in terms of pure gas is 
clearly unacceptable. This points easily to a further contribution 
from dust. According to our best fit, the total column density 
of gaseous iron is then JVpf = (0.56 + 0.08) x 10 16 cirT 2 . 
From the iron edge fit we obtain a total dust column density of 

Npe St = ( 3 - 8 ± °- 5 ) x 10 ' 6 cm_2 - ^ sum of ^ Gse iron compo- 
nents provides Np e = (4.3 +0.5) x 10 16 cm~ 2 . When we compare 
this number to the value predicted from our set of solar abun- 
dances, Afp™ = (5.10 + 0.06) x 10 16 cirT 2 (given the best fit total 
hydrogen column density of Nu = (1-63 + 0.02) x 10 21 cirT 2 ) 
we obtain that the abundance is 0.85 + 0.08 times solar. The 
depletion of iron, meant here as the ratio of dust over the total 
amount of a given element is instead 0.87 + 0.14 (Table |6}. 
This value is often reported to be higher (e.g. 0.97, [Jenkins 
2009). Howeve r, other stud ies report lower values for the Fe 
depletion (0.7, IWilms et al.l |2000|) . We note that, considering 
some uncertainties that still remain in the edge fitting (Sect. [43) 
the value we obtain may be considered as a lower limit. 
At a temperature of kT =0.5 eV, Fe is mainly Fei and only 
~2% of the total iron in Fen. The sum of the mildly ionized 
phase (Fe 11-Fe iv), excluding therefore the high ionization ions 
described in Sect. 14.11 is 0.24 x 10 16 crrT 2 along this line of 
sight, i.e. Fe; on /Fef a t s ~0.3, remembering that neutral iron is 
~ 87% depleted. 

In the oxygen region the evidence for dust is not as striking 
as in the iron region (Sect. 14.2b . From the best fit including 
cold (kT = 0.5 eV) gas, the amount of oxygen in the gas form 
is consistent to be solar (N^ s = (9.8 + 1.0) x 10 17 crrT 2 ). 
However, the amount of oxygen in dust compounds is 
jV ust = (2.5 + 0.2) x 10 17 crrT 2 . The total amount of oxygen is 
then No = (1.2 + 0.1) x 10 18 crrT 2 , showing about 23% over- 
abundance with respect to solar (Table [6j. With this reference 
value, the amount of depletion is 0.20 + 0.02. Oxygen depletion 
and abundan ces, based on Ch andra data have been previously 
presented bv lYao et al.l (120061) . As expected, the column density 
of Oi (Nq S ) that we measure is similar to their s. However, 
taking as a reference fhe lAnders & Grevessd(ll989l) abundances 
(A£ ?89 = 8.5 x lO^. lYao et ail (f2006b find an underabundance 
of oxygen of about 30% . In the reference list that we adopted 
(Lodders & Palmd [2009), the absolute oxygen abundance is 
significantly lower (A™ = 6.0 x 10~ 4 ). 

Finally, the amount of mildly ionized oxygen (On-Ov, with 
total column density 5.5 x 10 I6 crrT 2 ) over the total oxygen in 
gas form is Oi on /Of a t s ~0.05. 



As shown in Sect. 14.51 the inclusion of Si in dust is 85%, 
while for Mg is 97% of the total ISM. These values are com- 
parable to what has been reported in the literature (e.g. 80- 
92 and 90-97 % for Mg and Si respectively. IWilms et al.M2000t 
IWhittet ll2003l) . However, Mg and Si could not be studied in de- 
tail in this case because of the relatively low column density 
toward 4U 1820-30 which imprints shallow absorption edges. 
Moreover, our fitting of the Mg and Si edges relies on the oxygen 
modeling. Therefore, if other Mg or Si compounds are present, 
here they are difficult to detect. Keeping in mind these limita- 
tions we report the depletion values of Mg and Si in Table [6] for 
completeness. The abundance estimate of both Mg and Si are 
again driven by the oxygen edge fitting and are formally slightly 
above the solar values. 

In Table [6] we list for each element, the ratio with Proto-solar 
abundances and the amount of depletion defined as the ratio be- 
tween dust and total ISM abundance. 

A further test on the abundances derived from the X-ray data 
would of course come from high-resolution (R = A/AA > 
20,000) UV data, wh ere the elemental gas pha se can be ac- 
curately studied (e.g. I Savage & Sembachl 119961) . In the case 
of 4U 1820-30, high-resolution UV data (either from HST or 
FUSE) are not available. 4U 1820-30 has been observed only 
with the HST-STIS-G140L spectrograph (R = 1,000). We find 
that the signal to noise ratio of the data was insufficient for a 
quantitative study on the absorption lines. 

5.4. The location of the cold matter 

It has been established that a gradient in the abundances in our 
Galaxy exists for the most abundant metals. An average slope of 
0.06dexkpc~' should roughly apply to the gradient for both O 
and Fe (see Sect[T|). 

In our analysis we estimate a slight overabundance of oxygen (by 
a factor ~ 1.2, see Tabled, while iron is a factor ~ 0.85 of the 
solar value. This appears to be in contradiction with the expected 
abundances at the distance of the source, where according to 
the gradient above, some overabundance is expected. However, 
4U 1820-30 is located at latitude b = -7.9133°, i.e. about 1 kpc 
below the Galactic disk, therefore our line of sight intercepts a 
relatively small fraction of the cold ISM near the source. The 
cold phase that we detect towards this source is rather due to ab- 
sorption in the environment close to the Sun. 
For iron, a gradient as a function of both the radial distance 
from the Galac tic center and the height above the disk has been 
estimated (e.g. IChen et aT 2003, based on open clusters mea- 
surements). Our value of [Fe/H0 is -0.07. This is consistent, 
within the errors, with absorption in the disk at the distance of 
the Sun rather than absorption local to the sour ce, at the height 
of 4U 1820-30 below the disk (IChen et alj2003l) . 
The location of the gas based on the oxygen abundance is diffi- 
cult t o determine, as ther e is a large scatter in the measurements 
(e.g. lRudolph et alT2 006, for a compilation of results). Our oxy- 
gen abundance fits well with both absorption far or near the ob- 
server. It is therefore not straightforward to understand where 
the bulk of the absorption takes place on the basis of gradient 
measurements. However, a systematic analysis of dust scattering 
halos points out that most of the scattering (and therefore the ab- 
sorpti on by dust) should happen within ~ 3° from the Galactic 
plane (Predeh l & Schmittlll995l) . This further support the intu- 
itive idea that most of the cold absorbing material is located close 
to the observer. 



9 defined as log(Fe/H)-log(Fe/H), 
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Table 6. Relative abundances and depletion values from the 
present analysis. 



elem. 






it 


A z /A° 


dust/ISM 


N 


13.3 + 0.1 


0.0 




1.2 ±0.1 


0.0 





98 ± 10 


25: 


t2 


1.23 ± 0.03 


0.20 ± 0.02 


Fc 


0.56 ± 0.08 


3.8 


±0.5 


0.85 ± 0.08 


0.87 ±0.14 


Nc 


20.7 ± 0.2 


0.0 




0.9 ±0.1 


0.0 


Mg 


< 0.3 


7.8 


±0.8 


1.28 ±0.13 


>0.97 


Si 


< 1.3 


7.8 


±0.8 


i 9C+0.20 

1 - Z,J -0.12 


>0.86 



Notes: Abundances are referred to Lodders & Palme] i2009h 
Column densities are in units of 10 16 cm" 2 



5.5. The silicates 

In the present study, we find that the fit of both the iron and oxy- 
gen edges in general rejects silicate models containing iron (e.g. 
Fig- HIE]). In the iron region, the shift of the L III edge is con- 
sistent with absorption by metallic iron. Also a modest quantity 
of iron in the form of Fe2C>3 is allowed by the fit. As discussed 
above, in the oxygen region the evidence for dust is not striking. 
Our results may be still partially contaminated by instrumental 
effects, such as bad pixels in the RGS, which are however in- 
cluded in the response files using the most updated calibration. 
We keep limitations in mind when discussing the physical impli- 
cations of this result. 

In summary, we find that absorption by dust is mainly 
caused by metallic iron and glassy-enstatite. This kind of 
composition is reminiscent of the composition of GEMS (Glass 
with Embedded Metal and Sulfides), that are smal l grains 
abundant among the interplanetary dust particles (iBradlevi 
119941) . Most of GEMS particles should not have an ISM 
origin, but they rather res ide in interplanetary environments 
dKeller & Messengerfl2008h . However, GEMS with anomalous 
comp osition may have been processed in the ISM dMatzel et alj 
2008, and references therein). In particular, some of those 
particles show a low am ount of sulfur relative to silicon 
dKeller & Messenger] 120081 S/Si 0.19) , which i s more similar 
to what is found in the diffuse ISM (Sofia 2004) rather than in 
the solar neighborhood (Anders & Ebihara 1982). However, we 
cannot yet test the contribution of sulfides in a typical GEMS 
in our data, as FeS laboratory measurements around both the 
S K- and the Fe L-edge are not available. From the model, 
the glassy-enstatite MgSiO^ provides the best fit. Therefore 
the Mg/Si ratio is 1 by definition. This is also the value found 
for anom alous composition GEMS which are possibly of ISM 
origin, dKeller & Messenger 2008). Unfortunately we have not 
yet the means of testing compounds with a varying amount of 
Mg (or Fe). Therefore, we cannot test wether a Mg/Si ratio 
of ~0.6 (typical of the av e rage GEMS in in terplanetary dust 
iKeller & Messenger! l2004t llshii et alj l2008h would be still 
allowed by the fit. 

We find that the Fe/Si ratio ranges from 0.42-0.55, which 
matches this ratio in any type of G EMS where Fe/Si=0.43-0.54 
dKeller & Messengerfl2004T 120081) . This is differe nt from wha t 
is expected in the diffuse ISM (Fe/Si=0.8 + 0.2, ISofial I2004T) . 
However, the Fe/Si ratio is currently under debate. Recent 
studies have revised this v alue, lowering it in fact significantly 
to roughly Fe/Si~ 0.1 (see iMin et at]|2007l for a discussion). 
Finally, the amount of Fe along this line of sight is about two 
times less than Mg (Mg/Fe= 2.0 + 0.3). This estimate may be 
biassed because we cannot tune the amount of Mg, which is 
bound to be equal to Si in MgSi03. Another caveat that should 



be kept i n mind is that large grains would be grey to X-ray 
radiation dWhittet Il200"3t) . Therefore, at least in absorption, the 
large grain popula tion (which ma y contain more iron) can be 
under-represented dLee et al . 200§. 

Even taking into account the limitations imposed by our 
models, we find that, from both the O and Fe edges fits, Mg- 
rich rather than Fe-rich silicates are present along the line of 
sight to 4U 1820-30. Such a result stren gthens previous studies 
in both IR dMin et al.ll2007h and X-ravs (ICostantini et al]|2005l) 
which point to the same conclusion. Interestingly, a similar spec- 
tral analysis, but restricted to the Fe L edge only, shows that ox- 
ides, rather than olivines or pyroxene are r esponsible for th e ab- 
sorption along the line of sight to Cyg X- 1 dLee et al . 2009). The 
column density towards 4U 1820-30 is about 4.5 times smaller 
than for Cyg X-l. This result may therefore point out a chemical 
homogeneity on different path lengths within the diffuse ISM. 

5.6. A fast outflowing gas? 

The oxygen region of the RGS spectrum of 4U 1820-30 displays 
evidence, clearly detected in two set of RGS observations, of 
two absorption lines, which are consistent with O iv and O v out- 
flowing at v ~ 1200km s" 1 (Fig. |7]). In the 4U 1820-30 accretion 
flow heavy elements should not be abundant, as the companion 
is classified as a He- white dwarf, on the b asis of X-ray binary 
evolutionar y models dRappaport et al.lll987l) . From optical spec- 
tra analvsis liSfelemans et al .1 ( 120101) found that elements heavier 
than Ne are absent in ultracompact systems. Oxygen is a rela- 
tively light element which could still be present in the outer en- 
velope of the companion and be transferred in the accretion disk. 
Moreover, a way of producing oxygen could be through triple-a 
burning of He in the white dwarf. Qualitatively, this could jus- 
tify the detection of an outflow containing oxygen. However this 
result needs to be tested with additional observations. 
Here we fitted the absorption in terms of a photoionized gas. 
Other lines are predicted by the model, but they are too weak or 
in a noisy part of the spectrum to be significantly detected. This 
system is well detected in both RGS data sets, taken 8 years 
apart. The outflow did not change its physical parameters during 
that time. We note that this models could not be tested against the 
Chtmdra-HETG data, because of the noise affecting the oxygen 
region. The model predicts also a substantial amount of C iv. We 
qualitatively checked STIS-G140L low-resolution data (taken in 
1998, about 3.5 years before the first RGS measurement) for the 
presence of a blueshifted C iv doublet. Fixing the outflow veloc- 
ity to the X-ray value we obtain an upper limit for the C iv col- 
umn density of < 10 13 cirT 2 , which is almost two orders of mag- 
nitude lower than that required by the X-ray model. Therefore, 
in the hypothesis of an outflow, this must have been absent three 
and a half years before the first RGS pointing, when the STIS 
observation was carried out. Fast outflows are not unusual in 
X-ray binaries an d are interpreted as accretion disk winds (e.g. 
Miller et al. 2006). However, the ions involved in the outflow are 
generally more highly ionized (e.g. O vn, Ne ix-Ne x). 
In terms of goodness of fit, this outflowing system can be equally 
well fitted by a collisionally ionized plasma with temperature 
T = 13.9 ± 0.8 eV and N H = 3.4 + 0.5 x 10 19 cirT 2 . The amount 
of C iv predicted by this model would be ~ 6 x 10 14 cm -2 cm -2 , 
less than what is predicted by the photoionized gas, but still in 
disagreement with the upper limit derived from the STIS data. 
However, higher-quality STIS data are necessary for a reliable 
comparison between the UV and X-ray band outflowing ab- 
sorber. 
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Fast moving collisionally ionized clouds in the line of sight 
with such high vel ocities are not reported by the UV surveys 
dSavage et al.l l2004). Therefore the phenomenon should arise in 
the proximity of the source, possibly in a less ionized impact 
region, that would cause ionization by collisions. A candidate 
is the spot where the accretion flow impacts the ac cretion disk, 
which should be less ioni zed than the disk itself dBoirin et al.l 
l2005t Ivan Peet e t al. 2009). The absence of emission lines from 
this gas implies that, regardless of the absorbing mechanism, 
the flow pointing toward the observer must be very collimated. 
Although plausible, the scenarios depicted above need support 
from further (multiwavelength) observations. 
Finally, absorption by dust seems unlikely, as dust features are 
generally smoother (Appendix A). 

6. Conclusions 

In this paper we present the X-ray analysis of the continuum of 
4U 1820-30 (using a quasi-simultaneous observation of XMM- 
Newton and INTEGRAL) and of the absorption features due to 
the cold matter in the line of sight (using XMM-iVewfon-RGS 
and Chandra-MEG data). 

The continuum shape and the Eddington ratio show that 
the source was caught in a high-state. The continuum is well 
fitted by black body emission plus a Comptonization component 
which extends up to 40keV. We do not find evidence of iron 
emission, either from neutral or ionized matter. This may be 
naturally explained by the metal-poor accretion stream expected 
from the white-dwarf companion. 

The absorption spectrum shows the presence of many com- 
ponents with different ionization. We focused on the cold and 
mildly ionized phase only. Oxygen has been found slightly over- 
abundant by a factor 1 .23 times the solar value. Iron is on the 
contrary slightly underabundant (~0.85 times solar). The abun- 
dance values are not dramatically deviating from the solar ones 
and do not allow us to assign a precise location of the absorbing 
gas. However, a location close to the observer seems likely. 

Thanks to the simultaneous study of absorption by dust and 
gas we measured also the element depletion. Oxygen is mildly 
depleted by a factor about 0.20. The depletion of iron is more 
evident, as the depletion factor is 0.87. The depletion of Mg and 
Si are more difficult to determine. We find that they are depleted 
of a factor > 0.97 and > 0.86, respectively 

We modeled the dust contribution with the currently avail- 
ables absorption profiles of dust compounds. Our conclusions 
bear the uncertainty due to the still limited dust data-base and 
a lower sensitivity in selected spectral regions. However we 
clearly find that both the oxygen and iron edges cannot be 
fitted by iron-rich silicates. On the contrary, the oxygen edge is 
consistent to be mostly absorbed by enstatite (MgSi03, possibly 
in a glass-form). Metallic iron should be the main absorber in 
the iron edge. This leads to the interesting possibility that a 
GEMS-like form (Mg-rich silicates with metallic iron inclusion) 
of grain may be absorbing along this line of sight. A fraction 
of the studied GEMS, in particular the sulfur-poor grains, are 
believed to be of ISM origin and have also been proposed as 
constituents of ISM. For the first time an X-ray absorption 
analysis provides a tentative confirmation of this scenario. 

Finally, we report the tentative detection of a mildly ionized 
outflow (v out ~ 1200 km s~'), highlighted by the Oiv and Ov 



absorption lines. Both a photo- or collisional- ionizing process 
could fit the lines, leaving open the interpretation on the nature 
of this gas. 
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Appendix A: Oxygen compounds 

Here we show the profiles of the absorption around the oxy- 
gen edge for the compounds used in this analysis. These and 
other profiles are inclu ded in the AMOL model, implemented 
in SPEX. We refer to iPinto et al.l (1201 Ol) for a complete list 
and description of the oxygen compounds (with the exeption of 
MgSiOO. 



Fig. A.l. Transmission of the dust absorption models included 
in the present analysis for the oxygen region. The oxygen col- 
umn density has been se t here to 10 18 cm" 2 for all compounds. 
allBarrus et alJd 19971). bllParent etaT] (120021). cllvan Aken et all 
( 19981). d) iLee et alJ d2008h . See also fpinto et al) d2010h for de- 
tails. 
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